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Abstract 

The decays of J/ip — > ujKKtt and J/ip — > 4>KKtt are studied using 5.8 x 10 7 J/ip events collected 
with the Beijing Spectrometer (BESII) at the Beijing Electron- Positron Collider (BEPC). The 
KgK^-K^ and K + K~ir° systems, produced in J/ip — > ujKKtt, have enhancements in the invariant 
mass distributions at around 1.44 GeV/c 2 . However, there is no evidence for mass enhancements 
in the KKir system in J/ijj — > cpKKir. The branching fractions of J/ip — > uK^K^tt^ , cpKgK^ir^ , 
ujK*K + c.c., and cj)K*K + c.c. are obtained, and the J/ip — > nKgK ± 7r T branching fraction is 
measured for the first time. 

PACS numbers: 13.20.Gd, 13.25.Gv, 13.20.-v, 12.38.Qk, 14.40.-n 
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I. INTRODUCTION 



A pseudoscalar gluonium candidate, the so-called E/l( 1440), was observed in pp annihi- 
lation in 1967 [l| and in J/if> radiative decays in the 1980's [2,0,01. After 1990, more and 
more observations revealed the existence of two resonant structures around 1.44 GeV/c 2 in 
a (980)7r, KKti, and K*K spectra @, 0, 0, 0, 0, ffl, Ej. They showed that the lower state, 
77(1405), has large couplings to a (980)7r and KKti, while the high mass state, 77(1475), fa- 
vors K*K. The 77(1405) was also confirmed by Marklll P3, Crystal Barrel 0, and DM2 Q 
in decays into 777r7r in J / if> radiative decays and pp annihilations. 

In contrast, although 77(1475) was observed in KKn (K*K) @, 0, 0, 0, 0, 0, [H|, it was 



not seen in 77^71- . Nonetheless, the study of KKix and 777777 channels in 77 collisions [14 
showed that 77(1475) appeared in KKtc, but not in 777777, while r/(1405) appeared in neither 
channel. 

The study of the decays J /if) — > {^,uj,<^\KKt\ is a useful tool in the investigation of 
quark and possible gluonium content of the states around 1.44 GeV/c 2 . In this paper, 
we investigate the possible structure in the KKti final state in J /if) hadronic decays at 
around 1.44 GeV/c 2 , and measure the branching fraction of J /if} — > t]K s K ± tt t for the first 
time, based on 5.8 x 10 7 J/if> events collected with the Beijing Spectrometer at the Beijing 
Electron-Positron Collider (BEPC) . 



II. THE BES DETECTOR 



BESII is a large solid-angle magnetic spectrometer that is described in detail else- 
where [l5j. Charged particle momenta are determined with a resolution of u p /p = 
1.78%^/1+p 2 (p in GeV/c 2 ) in a 40- layer cylindrical main drift chamber (MDC). Par- 
ticle identification (PID) is accomplished using specific ionization (dE/dx) measurement 
in the MDC and time-of-flight (TOF) information in a barrel-like array of 48 scintillation 
counters. The dE/dx resolution is (JdE/dx — 8.0%; the TOF resolution is otof = 180 ps for 
Bhabha events. Outside of the time-of-flight counter is a 12-radiation-length barrel shower 
counter (BSC) comprised of gas proportional tubes interleaved with lead sheets. The BSC 
measures the energy and direction of photons with resolutions of &e/E ~ 21%y/E (E in 
GeV), o"0 = 7.9 mrad, and o z = 2.3 cm. The iron flux return of the magnet is instrumented 
with three double layers of counters that are used to identify mouns. 

A Geant3 based Monte Carlo (MC) package (SIMBES) with detailed consideration of 
the detector performance is used. The consistency between data and MC has been carefully 
checked in many high purity physics channels, and the agreement is reasonable [3]. The 
detection efficiencies and mass resolutions for each decay mode presented in this paper are 
obtained with uniform phase space MC generators. 



III. ANALYSIS 



In this analysis, 00 mesons are observed in the u — > 7r + 7r 77 decay, <f> mesons in the 
(f> — > K + K~ decay, and other mesons are detected in the decays: Kg — > 7r + 7r~, tt° — > 77, 
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7] — > 7T + 7r 7T°. The final states of the analyzed decays J /if) — > {uj,i]}KgK ± TT T , uoK + K 7r°, 
^K^K^ir*, and (f>K+K-7r° are 2(n + n-)K ± 7i T ^, tt+tt" fT + K+.?r-7r+7r-.Fr ± 7r T , 
and 2(K + K~) / ~f / y, respectively. 

Candidate events are required to satisfy the following common selection criteria: 



1. The correct number of charged tracks with net charge zero is required for each event. 
Each charged track should have a good helix fit in the MDC, and the polar angle 9 of 
each track in the MDC must satisfy | cos#| < 0.8. The event must originate from the 
collision point; tracks except ^ from Kg must satisfy \/ x 2 + y 2 < 2 cm, \z\ < 20 cm, 
where x, y, and z are the space coordinates of the point of closest approach of tracks 
to the beam axis. 

2. Candidate events should have at least the minimum number of isolated photons asso- 
ciated with the different final states. Each photon should have an energy deposit in 
the BSC greater than 50 MeV, the angle between the shower development direction 
and the photon emission direction less than 30°, and the angle between the photon 
and any charged track larger than 8°. 

3. For each charged track in an event, XpidC*) * s determined using both dE/dx and TOF 
information: 

X 2 pi D {i) = X 2 dE/dx(i) + XTOF(i), 

where % corresponds to the particle hypothesis. A charged track is identified as a K 
(it) if Xpid f° r the K (tt) hypothesis is less than those for the 7r and p (K and p) 
hypotheses. 

4. The selected events are subjected to four constraint kinematic fits (4C-fit), unless 
otherwise specified. When there are more than the minimum number of photons in an 
event, all combinations are tried, and the combination with the smallest \ 2 ls retained. 



The branching fraction is calculated using 

B (JI^^ X ) = - x N° b \B(X^YV (1) 

£j/*l>->X^Y X ^J/-4> x rS{A — > Y ) 

and the upper limit for a branching fraction is determined using 

B{JI%P ^ X)< ^x % xB^y)x(l-^)' (2) 



where, N b s is the number of events observed, N up is the upper limit on the number of the 



observed events at the 90% C.L. calculated using a Bayesian method [17j . e is the detection 
efficiency obtained from MC simulation, Nj/^ is the number of J /if} events, (5.77±0.27) x 10 7 
Tsl. a sys is the corresponding systematic error, and B(X —>■ Y) is the branching fraction, 



taken from the Particle Data Group (PDG) [17j], of the X intermediate state to the Y final 
state. 
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A. J /if) {u, rftKgK^ 



At least one charged track must be identified as a kaon using TOF and dE/dx information. 
If there is more than one kaon candidate, the assigned kaon is the one with the largest kaon 
weight. Candidate events are fitted kinematically using energy momentum conservation (4C- 
fit) under the 2(7r + 7r~)A ± 7r =F 77 hypothesis, and \ 2 < 25 is required. Each event is required 
to contain one K s meson with six possible 7r + 7r~ combinations to test for consistency with 
the Kg. Looping over all combinations, we select the one closest to the K$ mass, denoted 
as m n + n -, provided it is within 15 MeV/c 2 of the K$ mass. After K$ selection, the two 
remaining oppositely-charged pion combinations along with the two gammas are used to 
calculate m n + n -^. Figure CD (a) shows the scatter plot of m 77 versus m 7r + 7r - 77 with two 
possible entries per event, where clear r\ and uj signals are seen. 




FIG. 1: (a) The scatter plot of m 77 versus m^+^-^j, and (b) the 7r + 7r~77 invariant mass for 
J /if) — ► {u, ^K^K^tt^ 1 candidate events with two possible entries per event. The curves in (b) are 
the results of the fit described in the text, and the shaded histogram in (b) shows the normalized 
background estimated from the A^-sideband region (0.025 GeV/c 2 < \iti w + w - — m K o \ < 0.055 
GeV/c 2 ). 

The 7r + 7r~77 invariant mass distribution with two possible entries per event is shown 
in Fig. [T] (b), where rj and uj signals are apparent. The branching fractions of J /if) — > 
u>KgK ± 7c T and r]KgK ± 7r T are obtained by fitting this distribution. The backgrounds for 
J /if) — > uKgK ± TT T which contribute to the peak in the uj signal region mainly come from 
non-A<j events and events from J/if> — > uK^Kg that survive selection criteria. The number 
of background events from J /if) — > u K s K s is estimated from Monte-Carlo simulation to be 
less than 2 . Backgrounds for J /if) — » r]A^A ± 7r =F contributing to the peak in the r\ signal 
region mainly come from non-A° events and events from J/ if) decays into A* ^ (1430)° — > 
( K s tt°) ( K s f]) . The latter contribution is estimated to be less than one event from MC 
simulation. Non-A^ events from the A^ sideband region (0.025 GeV/c 2 < \m n + n - —^^01 < 
0.055 GeV/c 2 ) are shown in Fig. [1] (b) as the shaded histogram, the background events are 
19.2 ±15.6 uj and —4.1 ±7.0 r\ by fitting the distribution with possible signals and polynomial 
background. 

A fit to the m^+ n -jj distribution is performed by using u and rj Breit-Wigner (BW) 
functions folded with Gaussian resolution functions plus a quadratic polynomial, shown 
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as the curve in Fig. [TJ (b). The numbers of events in the u and 77 peaks are 1971.7 ± 
41.4 and 231.6 ± 23.1, respectively. Here, the background events in the decays of J/ip — > 
uK s K ir^ and J/ip — >• f]K s K ± ir T estimated above are not subtracted but are included in 
the background systematic error. The J/t/j ujK^K+tt* and J/ip rjKgK^ detection 
efficiencies are obtained from MC simulation to be 1.48% and 1.18%, respectively. The 
branching fractions are then determined as: 

B(J/ip -> ujKgK^) = (3.77 ± 0.08) x 1(T 3 , 

B{J/ip — > t]K s K ± tt t ) = (2.18 ± 0.22) x 10" 3 . 
Here the errors are statistical only. 



1. J/%b -> luK*K + c.c. -> luK^K^t 

To select the u signal, the mass combination with 7r + 7r~77 closest to the u mass is 
required to satisfy \m K +i r - T y — < 0.04 GeV/c 2 . Figure [2] shows the scatter plot of m K o^ 
versus rriKn for J/ip — > u;K°i^ ± 7r =F candidate events, where the events in the cross bands 
correspond to J/^> — > uK*K + c.c. 
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FIG. 2: The scatter plot of m K o^ versus m K ± n ^ for J/tp — > uoK^K^tt^ candidate events. 

Figure[3](a) shows the scatter plot of m,+ rT; versus m n + n -, and there is an accumulation 
of events in the u and K s cross bands. The combined mass spectrum of K s n T and K ± 7r Zf 
in the signal region (box 1 in Fig. 121(a)), which is defined as \m n + n - — m K o \ < 0.015 GeV/c 2 
and \m n + 1 ^- T1 — m u \ < 0.04 GeV/c 2 , is shown in Fig. [21 (b), where a clear K* signal is 
observed. The K* signal is fitted with a BW function folded with a Gaussian resolution 
function plus a third-order polynomial, and 1208.3 ± 93.3 K* events are obtained. 

Non-cu and uon-Kg backgrounds are studied using u and K s sideband events. Figure [21 (c) 
is the fitted Kit mass spectrum in the u sideband region Qm n + n - — m K o \ < 0.015 GeV/c 2 , 
0.06 GeV/c 2 < |m 7r + 7r - 77 — mj\ < 0.14 GeV/c 2 , shown as horizontal sideband boxes 2 
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in Fig. [3] (a)) and K° s sideband region (0.03 GeV/c 2 < \m K +- K - - m K o\ < 0.06 GeV/c 2 , 
\m n + n - iy — m u \ < 0.04 GeV/c 2 , shown as vertical sideband boxes 3), and the number of K* 
sideband events iV^i = (686.2±56.0) is obtained. Figure[3](d) is background from the corner 
region (0.03 GeV/c 2 < \m^+ n - - m K o \ < 0.06 GeV/c 2 , 0.06 GeV/c 2 < \m n + n -^ - m w \ < 
0.14 GeV/c 2 , shown as diagonal boxes 4), and the number of K* events N sid 2 is equal to 
(134.1±25.5). The number of background events in the signal region is half of the sum of K* 
events in the u sideband and Kg sideband regions (AT si(il ) minus a quarter of the K* events 
in the corner regions (N sid2 ). So N bg = (686.2 ± 56.0)/2 - (134.1 ± 25.5) /4 = (309.6 ± 28.8). 





m(K7r) (GeV/c 2 ) m(Krv) (GeV/c 2 ) 

FIG. 3: (a) The scatter plot of m 7r + w -^ versus m n + w -, and the combined mass spectrum of K^ir^ 1 
and K^ir^ with two entries per event J/ifj — > ujK*K + c.c. candidate events for (b) the signal 
region (the central box 1); (c) the u and sideband regions (two horizontal boxes 2 and two 
vertical sideband boxes 3); and for (d) the corner region (four diagonal boxes 4). The curves are 
the results of the fit described in the text. 

The detection efficiency is estimated to be 1.23% from MC simulation. After background 
subtraction, the branching fraction is determined to be 

B(J/ip -^vK*K + c.c.) = (6.20 ± 0.68) x 10~ 3 , 

where the error is statistical only. 
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2. J/tp -> loX(1U0) 



Figure H] (a) shows the scatter plot of ni K o K ± w ^ versus m n + n -^, and Fig. 0] (b) is the 
K s K ± tc t invariant mass spectrum after ui selection (|m 7r +, r - 77 — m u \ < 0.04 GeV/c 2 ). Figs. 
IU (a) and (b) show a resonance near 1.44 GeV/c 2 , denoted as X(1440). To ensure that this 
peak is not due to background, we have made studies of potential background processes 
using both data and MC simulations. Non-cu and non-iQj processes are studied with u and 
Kg mass sideband events, respectively. The main background channel J/if> — > u>2(7t + 7t~) 
and other background processes with 6-prong events are studied by MC simulation. In 
addition, we also checked for possible backgrounds with a ^MC sample of 60 x 10 6 J/ip — *> 



anything decays generated by the LUND-charm model [19|. None of these background 



sources produces a peak around 1.44 GeV/c 2 in the K S K ir^ invariant mass spectrum. 
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FIG. 4: (a) The scatter plot of m K o K ± 7T ^ versus m^^-yy and (b) the K s K^ic^ invariant mass 
distribution for J /if) — > loK s K^it^ candidate events. The curves in (b) are the results of the fit 
described in the text. 



The K s K ± ir^ invariant mass distribution is fitted with a BW function convoluted with a 
Gaussian mass resolution function (a = 7.44 MeV/c 2 ) to represent the X(1440) signal and 
a third-order polynomial background function. The mass and width obtained from the fit 
are M = 1437.6 ± 3.2 MeV/c 2 and T = 48.9 ± 9.0 MeV/c 2 , and the fit yields 248.8 ± 35.2 
events. 

Using the efficiency of 1.45% determined from a uniform phase space MC simulation, we 
obtain the branching fraction to be 



B(J/ip -> cjX(1440)) • B(X{UA0) -> KgKV) 
where the error is only the statistical error. 



(4.86 ±0.69) x 10" 



B. J ftp -> ujK+K-tt 



At least one charged track is required to be a kaon and the combined PID probability 
for K + K~7T + 7T~ is required to be greater than those for the K ± 7i T ii + 7i^ and 7r + 7r~7r + 7r~ 
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hypotheses. A 4C kinematic fit is made under the K + K^7T + 7i^A , y hypothesis. There are 
three combinations to form two 7r°'s, and further a six-constraint kinematic fit (6C-fit) with 
the smallest \\c * s ma de requiring two 7r°'s from four photons. Events with \\c < 50 and 
Xqc < 50 are selected. To reject the possible multiple photon background events, xlc * s 
required to be less than those for the K + K~n + 'K~2 r ) ) K + K~7T + 7T~3 , y, and K + K~-K + 7r~5 , y hy- 
potheses. Background events with K s decays, such as K* Kl (1430)° — > KgK ± n T {7i°, 2n }, 
and ~fK*K* — > / ~fK s K ± TT T TT°, are eliminated by requiring \m n + n - — m K o\ > 0.02 GeV/c 2 in 
the 7r + 7r~ invariant mass. 

There are two 7r + 7r~7r° mass combinations, and the one closest to the u mass, denoted 
as Tfi lx +- K — , is selected. The scatter plot of rriK+K-ir" versus m^+^-^o is shown in Fig. [5] 
(a), where the circle indicates some enhancement from J/ip — > cuX(1440) events in the 
LdK + K~7i° decay. 




FIG. 5: (a) The scatter plot of m^+^-- 7r o versus m 7T + n - 7T o, (b) the K^tt invariant mass distribution 
with two possible entries per event, and (c) the K + K~tt° invariant mass distribution for J/ip — > 
■k + it~it®K + K~it® candidate events. The curves are the results of the fit described in the text, and 
the shaded histogram (b) shows the normalized background estimated from the w-sideband region. 



1. J/ip -» ujK^KT ojK+K-tt 

To suppress the main K*° backgrounds, \m K ± 7r ^ — m K *o\ > 0.05 GeV/c 2 is required. In 
addition to the above selection, the further requirement of \m 7T +- K --K° — m u \ < 0.04 GeV/c 2 
is imposed. The combined mass spectrum of K + 7t° and K~ir° is shown in Fig. [5] (b), where 
the K*^ signal is seen, and is fitted to obtain the branching fraction of J/ip — > ujK^K^. 

Background events for ujK^K^ which could contribute to the peak in the K*^ signal 
region mainly come from events with K* decays, such as J/ip — > K^K^ (1430)° into 4-prong 
plus multiple photons, J/ip — » <pK*K, and J/ip — > jK*!^*, but their contributions can 
be ignored according to MC studies. It is further confirmed that the J/ip — > ujK^K^ 
background is negligible using uo and ir° sideband events. 

The K ±r K° invariant mass distribution in Fig. [5] (b) (2 entries/event) is fitted with a 
K*^ BW function with the mass and width fixed to PDG values [l7| plus a third-order 
polynomial. The number of K*^ events obtained is (175.6 ± 27.4). The detection efficiency 
is 0.32%, and the branching fraction of J/ip — > uK*K + c.c. is determined to be 
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B(J/ip ^uK*K + c.c.) = (6.53 ± 1.02) x 10" 3 , 
where the error is statistical only. 

2. J/ip -> luX(U40) -> loK+K-tt 

Figure [5] (c) shows the K + K~7T° invariant mass recoiling against the u, where a X(1440) 
signal is observed. We have also studied potential background processes using both data 
and MC simulations. Non-u; processes are studied with the u mass sideband events (0.06 
GeV/c 2 < Im^+Tr-TrO — m u \ < 0.10 GeV/c 2 ). Background with u decays is studied by MC 
simulations, similar to those of J/ip — > ujK*K + c.c. — > ujK + K~n°. In addition, we also 
checked for possible backgrounds using a MC sample of 60 x 10 6 J/ip — > anything decays 
generated by the LUND-charm model. In each case, the K + K~tt° mass distribution shows 
no evidence of an enhancement near 1440 MeV / c 2 . 

By fitting the K + K~tt° mass spectrum in Fig. [5] (c) with a B W function convoluted 
with a Gaussian mass resolution function (a = 14.2 MeV/c 2 ) plus a third-order polynomial 
background function, the mass and width of M = 1445.9 ± 5.7 MeV/c 2 and T = 34.2 ± 18.5 
MeV/c 2 are obtained, and the number of events from the fit is 62.1 ± 18.3. A fit without a 
BW signal function returns a value of —2 InL larger than the nominal fit by 31.7 with three 
degrees of freedom (d.o.f.), corresponding to a statistical significance of 5.0 a for the signal. 

The efficiency is determined to be 0.64% from a phase space MC simulation, and the 
branching fraction is 

B{J/ip -> cjX(1440)) ■ B(X(UA0) -> K+K'tt ) = (1.92 ± 0.57) x 10" 4 , 
where the error is statistical. 



C. J/ip -» ^K^K+it* 

Events with six charged tracks are selected, and at least two charged tracks must be 
identified as kaons. If there are more than two kaons, the two kaons with the largest kaon 
PID probabilities are regarded as the real kaons. The other charged tracks are assumed, 
one at a time, to be a kaon, while the other three to be pions, and these combinations 
of three kaons and three pions are kinematically fitted. The hypothesis with the smallest 
X 2 is considered as the right combination, and \ 2 < 20 is required. Two combinations of 
oppositely charged pions are used to reconstruct the signal, and the one closest to the 
Kg mass is required to be within 15 MeV/c 2 . 

The invariant mass of the two mass combinations formed with oppositely charged kaons 
are shown in Fig. [6J where a clear signal is observed. A fit to the K + K~ mass distri- 
bution in Fig [6] is performed to obtain the number of J/ip — > cj)K s K ± n T events. Back- 
grounds contributing to the signal peak mainly come from J/ip into 0/2(1525) — > <pnn, 
0r/ — > (pn-K + -K~ , (pK s K s , and 02(7r + 7r~) (excluding cpKgKg). From MC simulations of these 
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background channels, the number of background events in the signal region is less than 
one, and i^-sideband events also show that the background is negligible. 



"u 80 
> 




m(K + K") (GeV/c 2 ) 

FIG. 6: The K + K~ invariant mass distribution for J/t/i — ► K + K~ir + Tr~ K^ir^ candidate events 
with two possible entries per event. The curves are the results of the fit described in the text. 

The K + K~ mass distribution in Fig. [6] is fitted with a BW function convoluted with 
a Gaussian mass resolution function (a = 2.93 MeV/c 2 ) plus a third-order polynomial 
background function. The number of <ft events from the fit is 227.1 ± 19.0. Using the 
detection efficiency of 1.56%, the corresponding branching fraction is 

B(J/ij) -> <j>K° s K ± i<*) = (7.37 ±0.62) x 10" 4 , 
where the error is statistical only. 

1. J/ip -> (f>K*K + c.c. -> <t>K G s K±^ 

To remove most non-0 background events, the K + K~ combination closest to the (f) mass 
is required to satisfy \m K + K - — m<t\ < 0.015 GeV/c 2 . The scatter plot of m K o w versus 
for candidate events is shown in Fig. [7] (a), where the events in the cross band correspond 
to the K* signal. 

The scatter plot of m^+TT- versus ttik+k- is shown in Figure [7] (b), and there is an 
accumulation of events in the <fi and Kg cross bands. Figure [7J (c) shows the combined 
Kg7c T and K ± tt t mass spectrum for events in the signal region (box 1 in Fig. [7J (b)), which 
is defined as \m^+ n - — m K o\ < 0.015 GeV/c 2 and \m K + K - — m^\ < 0.015 GeV/c 2 . A fit 
yields 194.8 ± 25.0 K* events. 

The same method as used for the J/if> — > ujK*K + c.c. — > uoK^K ± t^ t analysis is used 
to estimate the number of background events in the signal region, and Nb g = (10.0 ± 6.6), 
which is neglected in the branching fraction determination. 

The detection efficiency of J ftp — > <fiK*K + c.c. in this decay is 1.42%, and its branching 
fraction is determined to be 

B(J/ip -»■ <pK*K + c.c.) = (2.08 ± 0.27) x 10" 3 , 
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FIG. 7: (a) The scatter plot of m K o n versus rriKn, (b) the scatter plot of m,+ r versus m K + K ~, 
and (c) the combined Kgir T and ir^vr^ invariant mass distributions for events in the signal region 
(box 1) for J/i/j — ► </>if<j.K' ± 7r = ' = candidate events. The curves are the results of the fit described in 
the text. 

where the error is statistical only. 



2. J/ip -> (pX(U40) -» <frKgK tt^ 

The distribution of i^^K ± 7r T invariant mass recoiling against the signal is shown in Fig. 
H] (a), and there is no evidence for X(1440). The upper limit on the number of the observed 



events at the 90% CL. is 8.1 17]. The likelihood distribution and the 90% CL. limit are 
shown in Fig. [8] (b). The likelihood values for the number of events are obtained by fitting 
the KgK^ir 11 distributions with a X(1440) signal plus a third-order polynomial background. 
Its mass and width are fixed to those in the decay J/ip — > u;X(1440) — > uK^K^-k^ . 




iti(KsK7t) (GeV/c 2 ) Number of Events 

FIG. 8: (a) The K ^^ invariant mass recoiling against the <p, and (b) the number of events of 
X(1440). The curve in (a) is a third order polynomial to describe the background, and the observed 
number of events at the 90% confidence level using a Bayesian method is indicated by the arrow 
in (b). 

The detection efficiency is 2.53%, and the upper limit on the branching fraction at the 
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90% C.L. is: 



B(J/i/) -> 0X(144O) -> <pK° s K + 7c- + c.c.) < 1.93 x 10" 



(3) 



D. J/V -» (j)K + K-TT C 



At least three charged tracks must be identified as kaons. A 4C-fit is applied under the 
hypothesis J/ip — > / -f-f2(K + K~), and x 2 < 16 is required. To reject possible background 
events from J/ip — > / y2(K + K^), the x 2 of the 4C fit for J/ip — > 772(i^ + i^~) is required to 
be less than that for the r j2(K + K~) hypothesis. There are four possible ways to combine 
the oppositely charged kaons in forming the 0, and the K + K~ combination closest to the <fi 
mass is chosen. Figured (a) shows the scatter plot of m 77 versus m K + K ~, and clear and 



7T 



signals are seen. 
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FIG. 9: (a) The scatter plot of to 77 versus rrix+K-i an d (b) the K^ir invariant mass distribution 
for events in the signal region (box 1) for J/ip — > r y"{2(K + K~) candidate events with two entries 
per event. The curves are the results of the fit described in the text. 



1. J/ip -> 4>K*±K* -> 0K+K-1T 



Figure [9] (b) shows the K ±r K° combined mass spectrum for events in the signal region (box 
1 in Fig. 0(a) ), which is defined as \m K + K - — m^\ < 0.015 GeV/c 2 and |m 77 — m n o\ < 0.04 
GeV/c 2 , and a clear K signal is seen. It is fitted with a BW, whose mass and width are 
fixed to those of K in the PDG, plus a third-order polynomial. The number of K events 
from the fit is 277.8 ± 27.7. The sidebands are used as before to estimate the number of the 
corresponding background events in the signal region, and the result is Ni g = (40.1 ± 10.1). 

After subtracting the above background and incorporating the efficiency of 1.71% from 
MC simulation, the branching fraction is determined to be 



B(J/ip -»■ <f)K*K + c.c. 
where the error is statistical only. 



(2.96 ±0.37) x 10" 
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2. J/tp -> <pX(lUO) -» (pK+K-ir 

The distribution of K + K~it° invariant mass recoiling against the is shown in Fig. [10] 
(a). No evidence for the X(1440) is observed near 1440 MeV/c 2 . The upper limit on the 
number of the observed events at the 90% CL. is 10.5 [l7j]. The likelihood distribution and 
the 90% CL. limit are shown in Fig. [TO] (b). The likelihood values for the number of events 
are obtained by fitting the K + K~tt° distributions with a X(1440) signal, whose mass and 
width are fixed to those of the decay J/if> — > ujK + K~tt° , plus a third-order background 
polynomial. The detection efficiency is 2.49%, and the upper limit on the branching fraction 
at the 90% CL. is: 

B(J/ip -> 0X(144O) -> (f)K + K-iT ) < 1.71 x 10~ 5 . (4) 




FIG. 10: (a) The K + K~tt° invariant mass recoiling against the (j), and (b) the number of events 
of X(1440). The curve in (a) is the third-order polynomial to describe the background, and the 
observed number of events at the 90% confidence level using a Bayesian method is indicated by 
arrow in (b). 



IV. SYSTEMATIC ERRORS 



In this analysis, the systematic errors on the branching fractions mainly come from fol- 
lowing sources: 

• MDC tracking efficiency 

The MDC tracking efficiency is measured from clean channels like J/ip —>■ A A and 
ip{2S) — > 7t + 7t~J/i/j with J/i/j — > fi + n~. It is found that the MC simulation agrees 
with data within 1% — 2% for each charged track. Therefore, 12% is taken as the 
systematic error on the tracking efficiency for the channels with six charged tracks 
and 8% for the channels with four charged tracks in the final states. 
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Photon detection efficiency 



The photon detection efficiency is studied from J/ip — > p°7r° events [20]. The results 
indicate that the difference between data and MC simulation is less than 2% for each 
photon. Therefore, 4% is taken to be the systematic error from the photon efficiency 
for the channels with two photons and 8% for the channels with four photons in the 
final states. 



PID 



The PID efficiency of the kaon is studied with J/ip — > K + K~tt° events. The average 
PID efficiency difference between data and MC is found to be less than 2%. In this 
paper, 2%, 4%, and 6% are conservatively taken as the systematic errors on PID 
efficiency for the channels with one, two, and three identified kaons, respectively. 



Kg reconstruction 

The K s secondary vertex reconstruction is checked using J/ijj — > A'* ± i^ =F (i ; r* ± — > 
iT^7r ± ) events. It is found that the difference of the efficiency between data and MC 
simulation is 2.8%, which is taken to be the systematic error from the K s secondary 
vertex reconstruction. 



Intermediate decay branching fractions 



The branching fractions for r\ 



from the PDG 17 



7T + 7T 7T°, U 



7T + 7T 7T°, 



and 



K + K are taken 



and the errors on these branching fractions are included as 
systematic errors in our measurements. The error on the K s — > tt + tt~ branching 



fraction is neglected in this analysis. 



Kinematic fit 



Kinematic fits are used to reduce backgrounds. Using the same method as in Ref . [21] , 
the decay modes J/ip — > 3(7r + 7r _ )7r°, J/ip — * 2(7r + 7r~)7r°, and J/ip — > 3(7r + 7r~) 
are studied [22J in order to estimate the corresponding systematic error. The 
kinematic fit efficiency differences between data and MC are 5.5%, 4.3%, and 8.7%, 
respectively. The efficiency difference between data and MC for the 6C-kinematic fit 



to rjc — > ojuj — > 2(n + ir n°) is about 10% [23j. Since the decays in this analysis are 



similar to the above decays, these systematic errors are taken as the corresponding 
systematic errors. 



Background uncertainty 
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TABLE I: Systematic errors in B{J/ip — ► {rj, u>, (f)}KKir). 



T lull 


TfiY c<Pl 71 


LOJ\ gJ\ 7T 


uK*K -> 


ujK*K -> 




A If* If 

(pj\ J\ — r 


<pl\ 1\ r 










ojK+K-it 






(JJXV XV /( 


H/I I Ui bU UI 






relative error ( 








A/TDP frnnVincr 
IVlL-'vy Li dCKlIl^ 


1 9 


1 9 


12 
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1 9 


1 9 


Q 
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JJilU LUil cllldtrilCy 
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A 


4 
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A 


JrdiLldc 11./ 


9 


9 


2 


2 


A 


A 


u 


i\ g znu. vei Lex 


2.8 


2.8 


2.8 




2.8 


2.8 




intermediate decays 


1.8 


0.8 


0.8 


0.8 


1.2 


1.2 


1.2 


kinematic fit 


5.5 


5.5 


5.5 


10 


8.2 


8.2 


4.3 


Back, uncertainty 


1.8 


2.0 


5.9 


11.7 


8.4 


6.1 


7.3 


MC statistic 


1.9 


1.3 


2.3 


3.2 


3.2 


3.4 


2.2 


MC model 






5.1 


5.1 




5.1 


5.1 


Number of J ftp events 


4.7 


total Sum 


15.4 


15.2 


17.1 


20.7 


18.4 


18.3 


15.6 



The background uncertainties come from the uncertainties associated with the estima- 
tion of the sideband backgrounds, the events from other background channels, as well 
as the uncertainties of background shapes, different fit ranges, and different binning. 
Therefore, the statistical error in the estimated number of background events, the 
largest difference from changing background shape, the difference from changing the 
fit range, the difference of changing fit binning, and some ignored backgrounds events 
are taken as the systematic errors due to the background uncertainty. 



• MC generator 



There may be interference between charged and neutral K* modes, which is not 
included in the MC generator. In the sample of J/ip — > ujK*K + c.c. — > uK^K^n^ 
decays, the K^n^ mass distribution in the K*^ sideband and signal regions and the 
KgTT mass distributions in the K*° sideband and signal regions from the scatter plot 
of mx±7r=F versus m^o,^, are studied in real data and MC simulation. It is found 
that the difference between data and MC sample is 5.1%, so 5.1% is taken as the 
systematic error from the MC model. 



Number of J/ip events 

The number of J/ip is (57.7 ± 2.7) x 10 6 , determined from J/ip inclusive four-prong 
events 18|. The uncertainty is taken as a systematic error in the branching ratio 



measurement. 



Tableland Table ILT1 list the systematic errors from all above sources, and the total systematic 
error is the sum of them added in quadrature. 
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TABLE II: Systematic errors in B{J/if) -» {u, (/>}A(1440) -> {w, 0}AAtt). 



J /ih 

j/yj — >■ 


wX(1440) 


uX(1440) 


0X(144O) 


c/>X(1440) 




-► uK^K^ 


-> L0K + K-TT° 




-> <j)K + K-ir° 


rjiioi bouice 


relative error (%) 




1 9 
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piioLon eiiiciency 
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2.8 




intermediate decays 


0.8 


0.8 


1.3 


1.2 


kinematic fit 


5.5 


10 


8.2 


4.3 


MC statistic 


2.7 


2.6 


0.8 


0.8 


Back, uncertainty 


6.4 


10.9 






Number of J /if) events 




4 


7 




Sum 


16.6 


19.6 


16.2 


12.6 



TABLE III: The branching fractions of J/^ decays in BESII. 



Decay 


final state 


No. of events 


efficiency 


Branching fraction (10 4 ) 


uA£A' + 7r~ + c.c. 


(vr + 7r-7r )AgA ± 7r T 


1971.7 ±41.4 


1.48% 


37.7 ±0.8 ±5.8 


7]K^K + ir~ + c.c. 


(Tr+vr-TT^AOA^T 


231.6 ±23.1 


1.18% 


21.8 ±2.2 ±3.4 


wA*A + c.c. 


(Tr+TT-TT^AOA^T 

(7r + 7r-7r°)A+A-7r° 


898.7 ±97.7 
175.6 ± 27.4 


1.23% 
0.32% 


62.0 ± 6.8 ± 10.6 
65.3 ± 10.2 ± 13.5 


(/>A° A+7T" + c.c. 


(A+A^AgA^ 


227.1 ± 19.0 


1.56% 


7.4 ± 0.6 ± 1.4 


4>K*K + c.c. 


(A+A-)A^A ± 7r =F 
(K+K )K + K'tt° 


194.8 ±25.0 
237.7 ±29.5 


1.42% 
1.71% 


20.8 ±2.7 ±3.9 
29.6 ±3.7 ±4.7 



V. RESULTS 



Table [ITT] lists the branching fractions of J /if) -> {rj, w, <f>}K%K ± it^ i J /if) -> {u, <p}K*K + 
c.c. from different decay modes. These branching fractions are somewhat larger than those 
of other experiments in Table [IV] [13] lEf but they are still consistent within errors. The 
branching fraction for J /if) — > T]KgK ir^ is measured for the first time. In the invariant mass 
spectra of KgK ± 7t T and K + K~tc° recoiling against the u, the resonance at 1.44 GeV/c 2 
is observed, with the mass, width, and branching fractions listed in Table 12 while in 
the invariant mass spectra of KgK ± ir T and K + K~n° recoiling against the <f>, no significant 
structure near 1.44 GeV/c 2 is seen and an upper limits on the J /if) decay branching fractions 
at the 90% C.L. are given in Table M 
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TABLE IV: The branching fractions of J/tp decays from Marklll Q and DM2 Q Collaborations 





Decay 


final state 


Branching fraction (10 4 ) 


Mar kill 


ujK^K + tt~ + c.c. 




29.5 ± 1.4 ± 7.0 




luK*K + c.c. 


(-k + -k--k°)K + K--k° 


53 ± 14 ± 14 




cj)K°K+-7r- + c.c. 


(K + K-)K° s K + t[- + c.c. 


7.0 ± 0.6 ± 1.0 


DM2 


(j)K s K + -K- + c.c. 


{K + K-)K^K ± it^ 


7.4 ± 0.9 ± 1.1 




4>K*K + c.c. 


\k + K-)K^K ± t^ 


20.8 ±2.7 ±3.7 



TABLE V: The mass, width, and branching fractions of J/ip decays into {u>, </>}X(1440). 



J/ip -> wA(1440) 

(X -> Ag A+Tr" ± c.c.) 


J/^f u;A(1440) 
(A -► AT+if-TT ) 


M = 1437.6 ± 3.2 MeV/c 2 
r = 48.9 ± 9.0 MeV/c 2 


M = 1445.9 ± 5.7 MeV/c 2 
r = 34.2 ± 18.5 MeV/c 2 



B(J/^> -* wA(1440) -* ujK q s K + tt- + c.c.) = (4.86 ± 0.69 ± 0.81) x 10" 4 
B(J/^ -> wA(1440) -> ujK + K--k g ) = (1.92 ± 0.57 ± 0.38) x 10~ 4 

B(J/j) ->■ </>A(1440) ->■ ^AgA+TT" + c.c.) < 1.93 x 10~ 5 (90% C.L.) 
B{J/i) -> <t>X(lte0) -> 0K+K-7T ) < 1.71 x 10" 5 (90% C.L.) 
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